THE DECOMPOSITION OF THE HYPERMETRIC CONE INTO 

L-DOMAINS 



MATHIEU DUTOUR SIKIRIC AND VIATCHESLAV GRISHUKHIN 



Abstract. The hypermctric cone HYP n+ i is the parameter space of basic 
Delaunay polytopcs of n-dimcnsional lattice. If one fixes one Delaunay poly- 
tope of the lattice then there are only a finite number of possibilities for the 
full Delaunay tessellations. So, the cone HYP„ + i is the union of a finite set of 
L-domains, i.e. of parameter space of full Delaunay tessellations. 

In this paper, we study this partition of the hypermetric cone into L- 
domains. In particular, we prove that the cone HYP n+ i of hypermetrics on 
n + 1 points contains exactly \n\ principal L-domains. We give a detailed 
description of the decomposition of HYP, l+1 for n = 2,3,4 and a computer 
result for n = 5 (see Tabic |2J). Remarkable properties of the root system D4 
are key for the decomposition of HYP5. 

1. Introduction 

An n- dimensional lattice L is a subgroup of ¥L n of the form L = v{L + ■ • • + v n 7j 
with («!,..., v n ) a basis of M n . Let S(c, r) be a sphere in R n with center c and 
radius r. Then, S(c, r) is said to be a Delaunay sphere in the lattice L if the 
following two conditions hold: 

(i) ||u — c|| > r for all v G L, 

(ii) the set S(c, r) D L has affine rank n + 1. 

The n-dimensional polytope P, which is defined as the convex hull of the set 
S(c,r) fl L, is called a Delaunay polytope of rank n. The Delaunay polytopes of 
rank n form a face-to-face tiling of W 1 . The Voronoi polytope Pv(L) of a lattice 
L is the set of points, whose closest element in L is 0. Its vertices are centers 
of Delaunay polytopes of L. The polytope Pv{L) forms a tiling of M. n under 
translation by L, i.e. it is a parallelohedron (see Figured]). 

The cones «S> , 5> are respectively the cone of positive definite, positive semi- 
definite n x n matrices. The rational closure S^ at>Q of iS> is defined as the 
positive semidefinite matrices, whose kernel is defined by rational equalities (see 
|DSV0 6]). Given a basis B = (i>i)i<i< n of a lattice L, we associate the Gram 
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The Delaunay polytopes of The Voronoi polytope of L 
L and its translates 

Figure 1. A lattice Lcl 2 and the induced partitions 

matrix a = B T B G 5™ . On the other hand if a G iS™ , then we can find an in- 
vertible real matrix B such that a = B T B, i.e. B is the basis of a lattice L = B7I 1 
with Gram matrix a. So, we can replace the study of Delaunay polytopes of L 
for the standard scalar product by the study of Delaunay polytopes of Z n for the 
scalar product x T ax. If one takes another basis B' of L, then B' = BP for some 
P G GL n (Z) and one has a' = P T aP, i.e. a and a' are arithmetically equivalent. 
In other words, the study of n-dimensional lattices up to isometric equivalence is 
the same as the study of positive definite n x n symmetric matrices, up to arith- 
metic equivalence. In [DSV06] it is proved that if a G S> , then one can define, 
possibly infinite, Delaunay polytopes of Z n for x T ax if and only if a G S^ at>0 . 

Given a polytope P of U 1 the condition that it is a Delaunay polytope for the 
norm x T ax translates to linear equalities and strict inequalities on the coefficients 
of a. An L-domain is the convex cone of all matrices a G iS" such that Z n has the 
same Delaunay tessellation for x T ax (Details see, for example, in |Vo08[ IDL97[ 
IDSV06] ). Voronoi proved that the cone S™ Q is partitioned into polyhedral L- 
domains. An L-domain of maximal dimension ^n(n + 1) is called primitive. An 
L-domain is primitive if and only if the Delaunay tiling related to it consists 
only of simplices. Each non-primitive L-domain is an open face of the closure 
of a primitive one. In particular, an extreme ray of the closure of an L-domain 
is a non-primitive one-dimensional L-domain. The group GL n (Z) acts on the 
L-domains of S™ by D t— > P T VP, and there is a finite number of orbits of L- 
domains, called L-types. The geometric viewpoint is most useful for thinking, 
and drawings about lattice and the Gram matrix viewpoint is the most suitable 
to machine computations. 

A metric on the set {0,1,..., n} is a function d such that d(x, x) = 0, d(x, y) = 
d(y, x) and d(x, y) < d(x, z) + d(z, y). A metric d is an hypermetric if it satisfies 
the inequalities 

(1) H z (d)= Yl z i z Ahj)<0 

0<i<j<n 

for all integral vectors z G Z n+1 such that Y17=o z % — ^-- The set of all hypermetrics 
on n points {0, . . . , n — 1} is denoted by HYP n . 
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n 


representative of orbits of facets of HYP n 


3 


(1,1,-1) (triangle inequality) 


4 


(1,1,-1,0) 


5 


(1, 1, —1, 0, 0) and (1, 1, 1, —1, —1) (pentagonal inequality) 


6 


(1,1,-1,0,0,0), (1,1,1,-1,-1,0), 




(1,1,1,1,-1,-2) and (2,1,1,-1,-1,-1) 


Table 1. The facets of HYP n for n < 6 



The group Sym(n) acts on HYP„; it is proved in |DGP06j that there is no 
other symmetries if n ^ 4. It is proved in [DL97j that HYP n+1 is polyhedral, 
i.e. among the infinite set of inequalities of the form ([1]), a finite number suffices 
to get all facets. This result can be proved in many different ways, see |DL97t 
Theorem 14.2.1]; the second proof uses that the image £(HYP n+ i) is the union of 
a finite number of L-domains. The purpose of this article is to investigate such 
decompositions of HYP n+1 . 

The set of orbits of facets of HYP n for n < 6 is given in Tabled] HYP 7 has 14 
orbits of facets (see [DL971 IBa99[ IRBj ) and the list is not known for n > 8. An 
inequality of ([T]) is called k-gonal if Y17=o \ Zi \ = ^- 3-gonal and 5-gonal inequalities 
are also called triangle and pentagonal inequalities, respectively. 

A Delaunay polytope P of a lattice L is called generating if the smallest, for 
the inclusion relation, lattice containing V(P) is L. Moreover, if there exist a 
family (vq, . . . ,v n ) of vertices of P such that for any v G L there exist ^ 6 Z 
with 

n n 
i=0 i=0 

then P is called basic and (vq, . . . ,v n ) is an affine basis. Given such an affine 
basis, we define the distance d(i,j) = \\vi — Vj\\ 2 and we have 

n 

H b (d)= b i b j d(i,j) = (r 2 -\\J2b i v i -c\\ 2 )<0, 

0<i<j<n i=0 

where r and c relate to the circumscribing sphere S(c, r) of P. So the hypermetric 
inequalities correspond to the inequalities determining a family (vq, . . . , v n ) to be 
an affine basis of a Delaunay polytope. Moreover we have Hf,(d) = if and only 
if Xir=o bi v i is a vertex of P. In other words the hypermetric cone HYP n+1 is the 
parameter space of a basic simplex in Z n . We refer for proofs to |As82[ IDL97j . 
In practice, if (vq, ■ ■ ■ , v n ) is an affine basis, we can replace it by vq = 0, Vi = e« 
and call the corresponding simplex main. At this point we should note that 
the hypermetric cone is just one possibility for a parameter space of Delaunay 
polytopes. Following |Er75llE~r92j . define M.2, n to be the space of real polynomials 
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of n variables with degree at most 2. We then have 

C n = {/ e M 2 , n | f(x) > for all x G Z n } 

If P is a Delaunay polytope of rank k < n, then we define 

C n {P) = {f eC n | f(x) = for all x G vert P} 

Those cones were used in |DR03j to find some so-called perfect Delaunay poly- 
topes. Note that if P = {0, e\, . . . , e n }, then the cone C n (P) is isomorphic to the 
cone HYP n+ i. 

The covariance map £ : d — > a transforms a hypermetric d on n + 1 points i, 
< z < n, of a set X into an n x n positive semidefinite symmetric matrix a as 
follows: 

<hi = Z(d(i,j)) = -(d(0,i) + d{0,j) - d{i,j)) 

(see |DL97l Section 5.2]). The covariance £ maps the hypermetric cone HYP n+ i 
into <S" at>0 . Note that there are n + 1 distinct such maps depending on which 
point of {0, ei, . . . , e n } is chosen as the zero point. 

2. Decomposition methods 

Recall that the Delaunay tiling related to a Gram matrix a from a primitive 
L-domain D consists of simplices. The set of Delaunay simplices of the tiling 
containing the common lattice point is the star Sto. By translations, along Z n , 
the star Sto determines fully the Delaunay tiling of Z n . The primitive L-domain 
T> belongs to £(HYP n+ i) if and only if its star Sto contains a main simplex. 

A wall W is an n( -" 2 +1 ' > — 1 dimensional L-domain, which necessarily separates 
two primitive L-domains L>, T>' . Let one moves a point a from the primitive L- 
domain T> to V by passing through W. When a G W, some pairs of simplices of 
Sto, which are mutually adjacent by a facet, glue into repartitioning polytopes. It 
is well known (see, for example, [DSV06t IVo08j ) that an n-dimensional polytope 
with n + 2 vertices can be triangulated in exactly two ways. When the point 
a goes from W into the L-domain T>' , each repartitioning polytope repartitions 
into its other set of simplices. 

Since each repartitioning polytope has n + 2 vertices, there is an affine depen- 
dence between its vertices. This affine dependence generates a linear equality 
between the coefficients of the Gram matrix a G «S" . This equality is just 
the equation determining the hyperplane supporting the wall W. If the point a 
lies inside the L-domain L>, then this equality holds as an inequality. 

So, the convex hull of vertices of any pair of simplices of Sto adjacent by a facet 
is a putative repartitioning polytope giving an inequality separating the L-domain 
T> from another L-domain. All adjacent pairs of simplices of Sto determine a 
system of inequalities describing the polyhedral cone of the primitive L-domain 
T>. Note that some of these inequalities define faces of T> but not walls. If the 
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n 


# primitive 


# facets 


# orbits of primitive 


# primitive L-domains 




L- types 


of HYP n+1 


L-domains in HYP n+1 


in HYP n+1 








under Sym(n + 1) 




2 


1 


3 


1 


1 


3 


1 


12 


1 


3 


4 


3 


40 


5 


172 


5 


222 


210 


8287 


5 338 650 



Table 2. Decomposition of HYP n+1 into L-domains 



adjacent pair of simplices contains the main simplex then the corresponding wall 
lies on a facet of the cone £(HYP n+ i). 

Using the above system of inequalities, one can define all extreme rays of T>, 
and then all facets of T>. The sum of Gram matrices lying on extreme rays of T> 
is an interior point a{T>) of D uniquely related to this L-domain (see |DSV06] for 
more details). Hence, L-domains T> and D' belong to the same L-type if and only 
if a(V) = P T a(T>)P for some P G GL n (Z), i.e. a(T>') and a(T>) are arithmetically 
equivalent. 

The algorithm for enumerating primitive L-domains in £(HYP n+1 ) works as 
follows. One takes a primitive Gram matrix a G £(HYP n+ i). There is a standard 
algorithm which, for a given a G <S> , constructs its simplicial Delaunay tiling. 
(For example, one can take a from a principal L-domain, described in following 
sections) . 

Using the star Sto of the Delaunay tiling, the algorithm, for each pair of ad- 
jacent simplices, determines the corresponding inequality. By the system of ob- 
tained inequalities, the algorithm finds all extreme rays of the domain T> of a, 
computes the interior central ray a(T>) and finds all facets of T>. The L-domain 
T> is put in the list £ of primitive L-domains in £(HYP n+1 ). 

Let F be a facet of T> which does not lie on a facet of £(HYP n+ i). For each 
repartitioning polytope related to the facet F of L>, the algorithm finds another 
partition into simplices. This gives the Delaunay tiling of the primitive L-domain 
V, which is neighboring to T> by the facet F. The algorithm finds all extreme rays 
of V, the ray a(V) and tests if it is arithmetically equivalent to a{T>) for some 
T> from C. If not, one puts T>' in C. The algorithm stops when all neighboring 
L-domains are equivalent to ones in L. This algorithm is very similar to the one 
in |SY04j for the decomposition of the metric cone into T-domains and belongs 
to the class of graph traversal algorithms. 

We give some details of the partition of HYP n+ i into primitive L-domains in 
Table H 

We now expose another enumeration method of the orbits of primitive L- 
domains in HYP n+ i. Consider a primitive L-domain T>. The group Stab(L>) = 
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{P G GL n (Z) : P T a(T>)P = a(T>)} is a finite group, which permutes the trans- 
lation classes of simplices of the Delaunay decomposition of D. It splits the trans- 
lation classes of simplices into different orbits. Let S be a basic simplex in L>; if 
one choose the coordinates such that vert S — {0, e±, . . . , e n } then one obtains an 
L-domain L>s, whose image by £ _1 is included in HYP n +i. A permutation of the 
vertex set of S induces a permutation in HYP n+ i as well. So, two cones ^(Vs) 
and ^(Vs') are equivalent under Sym(n + 1) if and only if S and S' belong 
to the same orbit of translation classes of simplices under Stab(P). Therefore 
from the list of L-types in dimension n, one obtains the orbits of L-domains in 
HYP n+1 . 



3. DlCINGS, RANK 1 EXTREME RAYS OF AN L-DOMAIN AND OF HYP n+ i 

We denote by b T c the scalar product of column vectors b and c. A vector 
v G Z n is called primitive if the greatest common divisor of its coefficients is 
1; such a vector defines a family of parallel hyperplanes v T x = a for a G Z. 
In the same way, a vector family V = (t>i)i<i<M of primitive vectors defines M 
families of parallel hyperplanes. A vector family V is called a lattice dicing if 
for any n independent vectors , . . . , Vi n G V the vertices of the hyperplane 
arrangement vj.x = cti form the lattice Z n (see an example on Figured]). This 
is equivalent to say that any n independent vectors V{ have determinant ±1, i.e. 
the vector family is unimodular. Given a dicing, the connected components of 
the complement of M n by the hyperplane arrangement form a partition of IR n by 
polytopes. It is proved in |ER94j that the polytopes of a lattice dicing defined by 
V = (t>j)i<j<M are Delaunay polytopes for the matrices belonging to the L-domain 
generated by the rank 1-forms (vivf)i<i<M, whose corresponding quadratic form 
is f Vi (x) = (vfx) 2 . The reverse is also proved there, i.e. any L-domain, whose 
extreme rays have rank 1 has its Delaunay tessellation being a lattice dicing. Such 
L-domains are called dicing domains; they are simplicial, i.e. their dimension is 
equal to their number of extreme rays. 

Our exposition of matroid theory is limited here to what is useful for the 
comprehension of the paper and we refer to [DG991 Aig79 ITr92] for more details. 



Given a graph G of vertex-set {1, . . . ,n} we associate to every edge e = 
a vector v e , which is equal to 1 in position i, —1 in position j and otherwise. 
The vector family V(G) = (v e ) ee E(G) is unimodular and is called the graphic 
unimodular system of the graph G. Given a unimodular n-dimensional system 
U of m vectors, for any basis B C U, we can write U = B(I n ,A), where A is 
a totally unimodular matrix and (I n , A) is the concatenation of I n and A. The 
matrix (— A T , I m -n) defines a unimodular system, which is called the dual of U 
and denoted Dual(C/). Given a graph G of vertex set V(G) and edge set E(G), 
we choose an orientation on every edge e and associate to it a coordinate x e . 
We define a vector space V to be the set of vector v G WL E ( G > satisfying for all 
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x G V(G) to the vertex cut equation 

with N(x) the neighbors of x and ei x ,y) = 1 if the orientation of the edge (x, y) 
goes from x to y and —1 otherwise. Take v±, . . . ,vn a basis of the space and 
denote by CoGr(G) the cographic unimodular system of the graph G defined to 
be the vector system obtained by taking the transpose of the matrix (vi, . . . , v^). 
The unimodular systems CoGr(G) and Dual(Gr(G)) are isomorphic. In |DG99] 
a general method for describing unimodular vector families is given using graphic, 
cographic unimodular systems and a special unimodular system named E§ (or 
Rio as in |Se80j ). 

Given a finite set E, a matroid M = M(E) is a family C(M) of subset of X 
called circuits such that: 

• for Ci, C 2 G C(M), it holds Ci £ C 2 , C 2 £ Ci if Ci ^ C 2 ; 

• if e G Ci n C 2 , then there is C 3 G C(M) such that C 3 C d U C 2 - {e}. 

A set of vectors g e , e G E, represents a matroid M(E) if, for any circuit C G 
C(M), the equality X] e ec ( ? e = ^ holds. A matroid, is called regular if it admits 
a representation as a unimodular system of vectors. If M(E) is a graphic or 
cographic matroid of a graph G with a set of edges, then circuits of M are 
cycles or cuts of G, respectively. 

A graph G is plane if it is embedded in the 2-plane such that any two edges are 
non-crossings. A plane graph defines a partition of the plane into faces delimited 
by edges. The dual graph G* is the graph defined by faces with an edge between 
two faces if they share an edge. Then (G*)* = G, and there is a bijection between 
(intersecting) edges of G and G* such that each cut of G corresponds to a cycle of 
G*, and vice versa. In other words, the cographic matroid of G and the graphic 
matroid of G* are isomorphic. 

The only rank 1 extreme rays of the cone HYP n+ i are cut metrics. For n < 5, 
the hypermetric cone HYP n+ i coincides with the cut cone CUT n+ i which is the 
cone hull of 2 n — 1 cut metrics. Denote iV = {1, . . . , n}; if S C N, S , ^0 then 
the cut metric 5s on X = {0} U iV is defined as follows: 

$s(hj) = 1 if H S\ = 1, and Ss(i,j) = 0, otherwise. 

The covariance map £ transforms the cut metric 8s into the following correlation 
matrix p(S) of rank 1: 

Pij{S) = 1 if {ij} C S, and Pij(S) = 0, otherwise, where 1 < i, j < n. 

The quadratic form corresponding to the correlation matrix p(S) is 

fs(x) = €(6 3 )(x) = ^p i:j (S)x i x j = C^Xi) 2 . 
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So, fs = f q with f g (x) = (q T x) 2 and q = J2 i&s h := b(S) the incidence vector of 
the set S. In summary: 

Lemma 1. A vector q G Z n determines an extreme ray f q o/£(HYP n+1 ) if and 
only if q = b(S) for some S'CiV, S ^ 0. 

By Lemma HJ if U determines a dicing domain in £(HYP n ), then the set of 
coordinates of vectors from U in the basis B forms a unimodular matrix with 
(0, l)-coemcients. Note that the columns of any (0, l)-matrix are incidence vec- 
tors of subsets of a set. Since any unimodular system U determines a dicing 
L-domain T>(U), we have the following proposition: 

Proposition 1. Let T>(U) be a dicing domain determined by a unimodular set 
U . The following assertions are equivalent: 

(i) V(U) lies in £(HYP n ); 

(ii) U is represented by a (0, 1) -matrix. 

4. The principal L-domain 

There is a unique primitive L-type, whose L-domains are dicing domains of 
maximal dimension \n{n + 1) ( [KoZo77l IDi72j ). Voronoi calls this L-type prin- 
cipal. Each principal L-domain is simplicial and all its \n(n + 1) extreme rays 
have rank 1. The set Q of vectors q determining extreme rays f q of a principal 
L-domain forms a maximal unimodular system. This system is the classical uni- 
modular root system A n representing the graphic matroid of the complete graph 
K n+1 on n + 1 vertices. 

In our case, when a principal L-domain is contained in £(HYP n+1 ), its extreme 
rays belong to the set {p(S) : S C N, S ^ 0} of extreme rays of the cone 
£(HYP n+ i). Hence, the vectors q have the form b(S) for S C N. We shall find 
all subsets of these vectors representing the graphic matroid of K n+ i. We orient 
edges of K n+ i into arcs and relate a vector b(S) to each arc of the directed graph 
K n+ i such that, for any directed circuit C in K n+1 , the following equality holds 

(2) X>&(S e ) = 0. 

Here b(S e ) is the vector related to the arc e and e e = 1 if the directions of e and 
C coincide, and e e = — 1, otherwise. 

Given a chain of equally directed arcs labeled by one-element set, a subchain 
of this chain determines a (0, 1 ^characteristic vector. It is known that the set 
of characteristic vectors of a set of connected subchains determines a graphic 
unimodular system. We show below that such graphical systems are contained 
in the set {b(S) : S C iV}. 

This relation of vectors b(S) and arcs of K n+ i provides a labeling of arcs of 
K n+ i by subsets S C N. We call this labeling feasible if the corresponding set 
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of vectors b(S) gives a representation of the graphical matroid of K n+1 , i.e. (j2J) 
holds for each circuit C of K n+1 . 

Consider a /c-circuit C = {e^ : 1 < i < k}, whose arcs have the same directions. 
Suppose that, for 1 < i < k, Si is a label of the arc e», and that this labeling is 
feasible. Then the equality Ya=i KSi) = holds. Since the coordinates of the 
vectors b(S) take (0, l)-values, this equality is not possible for a feasible labeling. 
Hence, the directed graph K n+1 with a feasible labeling has no circuit, whose arcs 
have the same directions. Any finite directed graph with no circuit has at least 
one source vertex (and a sink vertex as well). 

Now consider a directed 3-circuit C = {ei,e2,e 3 } of a feasible labeled K n+ i. 
Then two arcs of C, say the arcs ei,e2, have directions coinciding with the di- 
rection of C, and the third arc e 3 has opposite direction. If Si is a label of e^, 
i = 1,2,3, then we have the equality b(S\) + b(S 2 ) = b(S 3 ). This equality is 
possible only if Si fl S 2 = and Si U S 2 = S 3 . Since any two adjacent arcs of a 
complete graph belong to a 3-circuit, we obtain the following result: 

Lemma 2. Let two arcs and ej be adjacent in a feasible labeled graph K n+1 and 
have labels Si and Sj. Then Si fl Sj — if the directions of these arcs coincide in 
the 2-path [e*,ej]. If these directions are opposite, then either Si C Sj or Sj C Si. 

Let v be a source vertex of K n+ i and E(y) be the set of n arcs incident to v. 
Since all arcs of E{v) go out from v, any two arcs e, e' G ) have opposite 
directions in their 2-path [e,e'\. Let S(v) = {S e : e G -E'(f)}. By Lemma [2j the 
family is a nested family of n mutually embedded distinct subsets. This 
implies that the sets S G S(v) and the arcs e G i?(t>) can be indexed as Si, e^, 
1 < i < n, such that Si is the label of e, and = i. 

For 2 < z < n, let be the arc of the graph K n+1 , which forms a 3-circuit with 
the arcs ej_i,ej G E(v). Lemma [2] implies that the arc ^ has the one-element 
set Si — Si-i as label, and the direction of gi coincides with the direction of ej_i 
in their 2-path [ej_i,<7i]. Now, it is clear that the n arcs gi for 1 < i < n, where 
gi = ei, form an n-path, whose arcs have the same directions and are labeled by 
one-element sets. Recall that a non self-intersecting n-path in a graph with n + 1 
vertices is called a Hamiltonian path. We obtain the following result. 

Lemma 3. A feasible labeled complete directed graph K n+ i has a Hamiltonian 
path such that all its arcs have the same directions and each arc has a one-element 
labeling set. 

Let {0} U N be the set of vertices of K n+1 , where the vertex is the source. 
Let 0, ii, t2, ■ ■ ■ , i n be the vertices of the Hamiltonian path n in Lemma [21 The 
path 7r defines uniquely an orientation and a feasible labeling of K n+1 as follows. 
The arc with end- vertices ij,ik, where < j < k < n, is labeled by the set 
Sjk = {ir '■ 3 ' + 1 < r < k} C N. If one reverse the above order, then one gets the 
same family of sets, and the labeled graph K n+ i gives the same representation of 
the unimodular system A n . We have 
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Lemma 4. Any representation of the graphic matroid of the complete graph K n+1 
by vectors b(S), S C N, S 7^ 0, is determined by a complete order of the set N. 
Two opposite orders determine the same representation. 

Since there are n\ complete orders on an n-set, as a corollary of Lemma HJ we 
obtain our main result. 

Proposition 2. The cone £(HYP„ + i) contains \n\ distinct principal L- domains. 

So, each principal domain is determined by an order (and its reverse) of the 
set N. For the sake of definition, we choose the lexicographically minimal order 
O from these two orders. Let S{0) be the family of sets S C N, S ^ 0, such 
that elements of each set S determine a continuous subchain of the n-chain, 
corresponding to the order O. A principal domain determined by an order O of 
the set iV has \n(n + 1) extreme rays p(S) for S G S(0). 

Each face F of a dicing L-domain in £(HYP n+1 ) is uniquely determined by its 
extreme rays p(S), all of rank 1. Set 

S(F) = {S 1 C ]V : p{S) is an extreme ray of F}. 

Proposition 3. If n > A, then any two principal domains in £(HYP n+ i) are not 
contiguous by a facet. 

Proof. If two principal domains T>(0), V{0') C £(HYP n+1 ) share a facet F, then 
they have \n(n + 1) — 1 common extreme rays p(S) for S G <S(F). This implies 
that the families S(0) and S(0') should differ by one element only. But, for any 
two distinct orders O and O', the families S{0) and S{0') differ at least by two 
sets, if n > 4. For example, suppose O and O' differ by a transposition of two 
elements i and j. Then there is at least one subchain in O containing i and not 
containing j, which is not a subchain of O' . The same assertion is true for the 
order O'. □ 

5. The decompositions of HYP 3 and HYP 4 

Since, for n = 2 and n = 3 there exists only one primitive L-type, namely, the 
principal L-type, Proposition [2] describes completely the decompositions of the 
cones £(HYP n+ i) for n — 2,3. 

Recall that each facet F of HYP n+ i is described by an inequality (TjQ). If n < 4 
and F is a facet of HYP n+1 , then Zj e {0,±1}. Hence, we can denote triangle 
and pentagonal facets as F(ij\ k) and F(ijk; Im), respectively. Here Z{ = Zj = 1, 
Zk = —1, Zi — 0, I G X — {ijk}, for the triangle facet, and Z{ = Zj — — 1, 
Zi = z m = —1, z T = if r G X — {ijklm}, for the pentagonal facet. 

Note that £(<$s) = _p(S0- For 5 = {r/...A;}, set 5 = ij ...k and p(S') = 
p(ij ...A;). 

n = 2. The cone ^(HYPs) is three-dimensional and simplicial. There is only 
one order O = (12) with S(0) = {1, 2, 12}. Its three extreme rays p(l),p(2) and 
p(12) span a principal domain, i.e. £(HYPa) coincides with a principal domain. 
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n = 3. The cone £(HYP 4 ) is six-dimensional and has seven extreme rays p(i), 
i = 1,2,3, p(ij), ij = 12,13,23, p(123) and 12 facets £(F(ij;k)) for i,j,k G 
{0} U N. Note that 



Hence, the four-dimensional cone & = R+p(l) + K+p(2) + M+p(3) + M+p(123) 
intersects by the ray M. + p the three-dimensional cone C 2 = M. + p(12) + K + p(13) + 
IR + p(23). The ray M. + p is an interior ray of both. 

By Proposition (2j £(HYP 4 ) contains three principal domains. These three six- 
dimensional L-domains are determined by the three orders (123), (132) and (213) 
of N = {123}. Denote the domain determined by the order (ijk) as T>j, where 
j is the middle element of the order (ijk). Since the three one-element subsets 
{i}, i G N, and the set N give continuous chains in all the three orders, the four 
rays p(l) , p(2) , p(3) and p(123) are common rays of all the three domains T>i, 
T>2 and V 3 . Hence, the cone C\ is the common four-dimensional face of these 
three principal domains. The domain 2?j is the cone hull of C\ and two rays p(ij) 
and p(ik). The four triangle facets £(F(jk] i)),£(F(jk; 0)), £(F(0i; j)),£(F(0i; k)) 
of £(HYP 4 ) are also facets of ZV The other two facets of T>i separating the 
domain from Vj and V k are the cone hulls of C\ with the rays p(ij) and p(ik), 
respectively. 



A parallelohedron is an n-dimensional polytope, whose image under a transla- 
tion group forms a tiling of M n . Given a face F of a parallelohedron P, the set of 
faces of P which are translates of F is called the zone of P. For a parallelohedron 
P the Minkowski sum P + z(q) may not be a parallelohedron. A parallelohedron 
P is called free along a vector q and the vector q is called free for a parallelohedron 
P if the sum P + z(q) is a parallelohedron (see |Gr06j for more details on this 
notion) . 

If P q = P + z(q) is an n-dimensional parallelohedron, then P q has a non-zero 
width along the line l(q) spanned by q. This means that the intersection of P q 
with a line parallel to l(q) is distinct from a point. In this case, the lattice L q 
of the parallelohedron P q has a lamina H, i.e. a hyperplane H such that H is 
transversal to l(q), the intersection L q n H is an (n — l)-dimensional sub lattice 
of L q and each Delaunay polytope of L q lies between two neighbouring layers of 
L q parallel to L q PI H (see |DeG02j ). If a Voronoi polytope has a non-zero with 
along a line /, then the lamina H is orthogonal to /. 

6.1. Root lattice D 4 . The lattice D n is defined as 



p(l) + p(2) + p(3) + p(123) = p(U) + p(13) + p(23) = p. 



6. L-DOMAINS IN £(HYP 5 ) 



n 




n 




(mod 2)} 



i=i 
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If {ej : 1 < i < n} is an orthonormal basis of Z n , then the set of shortest vectors 
of D n is ±6j ± ej for 1 < i < j < n. It is the set of all facet vectors of Py{D n ) and 
form an irreducible root system, which we also denote by D n . There are three 
translation classes of Delaunay polytopes in D n : the cross polytope j3 n whose 
vertex set is formed by all e\ ± for 1 < i < n, the half cube \H n whose vertex 
set is {x G {0, l} n Xi = ^ (mod 2)} and a second half cube \H' n whose 

vertex set is {x e {1,2} x {0, l}™" 1 | J27=i x i = ( mod 2 )}- The two half 
cubes are equivalent under the automorphism group Aut(D n ) of the root lattice 
D n . It is proved in [Gr04t ]DcG03j that the Voronoi polytope Py(D n ) of the n- 
dimensional root lattice D n is free only along vectors which are parallel to edges 
of Py(D n ). Thus there are 2" + 2n free vectors. 

It turns out that when n = 4 all Delaunay polytopes are isometric to the 
cross-polytope /3 4 and equivalent under Aut(D 4 ). Any 2-face of /3 4 is contained 
in three Delaunay polytopes /3 4 of the Delaunay tessellation. This proves that 
the L-domains of D 4 are not dicing domains. Furthermore, D 4 is a rigid lattice 
(see |BaGr0lj ). i.e. its Delaunay tessellation determines the Gram matrix up to a 
scalar multiple. This means that D 4 determine a 1-dimensional L-type and that 
any L-domain containing it as an extreme ray is not a dicing domain. 

The free vectors of P v (D i ) are parallel and can be identified to the diagonals 
of the cross polytopes. They are 

• ±2e { with 1 < i < 4 for A 

• (±1, ±1, ±1, ±1) with even plus signs for |if 4 = /3 4 

• (±1, ±1, ±1, ±1) with odd plus signs for ^H' 4 = (3 4 

Up to a factor v^2, those 24 vectors are an isometric copy of the root system D 4 , 
which we denote by D 4 2 . The union D 4 U D 4 2 is the irreducible root system F 4 
(see [Hu90] l 

The Voronoi polytope Py(D 4 ) of the root lattice D 4 is the regular polytope 24- 
cell, whose automorphism group is the Coxeter group VK(F 4 ) of Schlafli symbol 
{3,4,3} (see |Cox48j ). Its number of vertices, 2-faces, 3-faces is 24, 96, 24. Each 
facet is an octahedron with four pairs of opposite and mutually parallel triangular 
2-faces. The facet vectors of the Voronoi polytope JV(D 4 ) are the 24 roots of the 
root system D 4 . The polytope Fy(D 4 ) has 12 edge zones of mutually parallel 
edges representing D 4 2 and 16 face zones of mutually parallel triangular faces. 
Each edge zone contains 8 parallel edges, and each face zone contains six parallel 
faces. 

We choose a basis of D 4 and denote by a(D 4 ) the Gram matrix of D 4 in this 
basis. The 24 vertices of Py(D 4 ) are given by |D 4) 2- 

6.2. The L-domains containing a(D 4 ). For n = 4, there are three primitive 
L-types of four- dimensional lattices: the principal type, and L-types called by 
Delaunay in [De29j types II and III. The ten-dimensional L-domains of these 
L-types are constructed as follows (cf., [DeG03] ). 
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Each /c-dimensional face of a principal domain relates to the graphic matroid of 
a subgraph G C K 5 on k edges. Hence each facet (of dimension 9) of a principal 
domain relates to the graphic matroid of K5 — 1, i.e. the complete graph K5 
without one edge. The Gram matrix a(D4) of the root lattice D4 is an extreme 
ray of L-domains of type II and III. The cone hull of a facet of a principal domain 
and of a ray of type a(D 4 ) is an L-domain of type II. Hence, any principal domain 
is contiguous in <S> by facets only with L-domains of type II. 

An L-domain of type II has the following three types of facets. One dicing 
facet by which it is contiguous in iS> to a principal domain relates to the graphic 
matroid of the graph K 5 — 1. Each of two other types of facets is the cone hull 
of the ray of type a(D 4 ) and a dicing 8-dimensional face related to the graphic 
matroids of K 5 — 2 x 1 or of K 5 — 2. Here each of these graphs is the complete 
graph K 5 without two non-adjacent or two adjacent edges, respectively. 

The complete bipartite graph K,., is formed by two blocks S\, S 2 of vertices with 
I Si I — i, \S 2 \ = j and two vertices adjacent if and only if they belong to different 
blocks. An L-domain of type III is the cone hull of a(D4) and a 9-dimensional 
dicing facet related to the cographic matroid CoGr( K33) of the bipartite graph 
K33. Each 8-element sub matroid of CoGr(K 33 ) is graphic and relates to the graph 
K 5 — 2 x 1. Hence, each other facet of an L-domain of type III is the cone hull of 
a(D 4 ) and a dicing 8-dimensional face related to K 5 — 2 x 1. In <S> , this facet is 
a common facet of L-domains of types II and III. 

So, an L-domain of type II is contiguous in iS> to L-domains of all three 
types. It is useful to note ( [RB051 IDV05] ) that if / belongs to the closure V of 
an L-domain of type II or III then the Voronoi polytope Pv(f) of Z n under the 
quadratic form / is an affine image of the Minkovski sum 

V(?) + AP y (D 4 ), A, > 0, A > 0, 

qeu 

where U is the unimodular set of vectors related to rank 1 extreme rays of T>, 
and Py(D 4 ) is the Voronoi polytope of the root lattice D 4 , whose form a(D 4 ) lies 
also on extreme ray of T>. 

6.3. Unimodular systems in D 4j 2- Let D(4) C D 4 2 be a subset of 12 roots 
chosen by one from each pair of opposite roots. The vector system D(4) is 
partitioned into three disjoint quadruples Qi, i = 1,2,3, of mutually orthogonal 
roots, given in Subsection 16. 1[ 

For what follows, we have to consider triples {ri,r 2 ,r 3 ) of roots chosen by one 
from each quadruple Qi, i.e. r { £ Qi, i = 1,2,3. Let t = (ri,r 2 ,r 3 ) be such a 
triple. Let {ijk} = {123}, i.e. these three indices are distinct. We have rf = 2, 
r I r j £ {^l}- Hence, the vectors 

(3) Tij =Ti - (rfr^rj for ij = 12, 23, 31, 

are roots of D 4) 2- Since r^r; = —rf^rj = 1, one of two opposite roots ±r»j belongs 
to Qk, say £ Qk- Hence r i2 £ Q3, r 2 3 £ Qi, r 3 i £ Q 2 . We have two cases: 
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(i) Tij belongs, up to sign, to the triple t, i.e. = r k , for all pairs ij; 

(ii) Tij does not belong to t, i.e. r^ ^ r k . 

In case (i), the vectors linearly dependent, and the triple t 

spans a 2- dimensional plane. We say that the triple t is of rank 2. Note that 
any two roots r, r' from distinct quadruples determine uniquely the third root 
r" = r — (r T r')r' and that there are 16 distinct triples of rank 2. 

In case (ii), the roots r, of the triple t are linearly independent. We say that 
the triple t has rank 3. In this case the roots rij for ij = 12,23,31, are distinct 
and do not coincide with the roots r,, i = 1,2,3. Moreover, it is not difficult to 
verify that the triple (ri2, r 2 3, r 3 i) has rank 2. 

Triples of rank 2 and 3 are realized in v^2-Py(D4) as follows. The roots of a 
triple t of rank 2 are parallel to edges of a 2-face of \ /f 2Py(D 4: ). We say that a 
triple t of rank 2 forms a face of v^Py(D4)- Since a triple of rank 2 forms a face 
of ^/2Py(D i ), the 16 triples of rank 2 relate to the 16 zones of triangular faces of 
it. 

Let t' = (?^,r 2 ,r 3 ) be a triple of rank 3. The 6 vectors ±r^, i = 1,2,3, have 
end-vertices in vertices of \f2P v {^>A). From each pair ±r ? - of opposite roots, we 
choose a vector such that rfrj = 1 for ij = 12, 23, 31. Then end- vertices of the 
roots Ti are vertices of a face F of \Z2Py(D 4 ), and the triple t = (r^ = r« — rj : 
ij = 12, 23, 31) has rank 2 and forms the face F. 

For what follows, we need subsets U of D(4) which are maximal by inclusion 
such that Py(D 4 ) + ^„ gC/ \z(q) is a parallelohedron. Note that each quadruple 
Qi is a basis of 1R 4 . Obviously, it is a unimodular set. It is a maximal by 
inclusion unimodular subset of D(A), since any other vector of D(4) has half- 
integer coordinates in this basis. However, it is proved in |DeG03| . that Py(D 4 ) + 
^2 q£U \z(q) is not a parallelohedron if U C D{4) is a quadruple, and it is a 
parallelohedron for any other maximal unimodular subsets U C D(4). Of course, 
a maximal unimodular set U ^ Qi does not contain each quadruple Qi as a 
subset. 

We show below that maximal unimodular subsets in D(4) represent either the 
graphic matroid of the graph K 5 — 1 or the cographic matroid of the graph K 33 
(see Figure [2]). 

The graph K5 — 1 is planar. Hence, the graphic matroid of K5 — 1 is isomorphic 
to the cographic matroid of the graph (K 5 — 1)*. Both graphs (K 5 — 1)* and K 33 
are cubic graphs on six vertices and nine edges. These graphs have a Hamiltonian 
6-circuit C$ on six vertices Vi, 1 < i < 6. Call Cq with the six edges (fj, fj+i) by 
a rim, and the other three edges by spokes. The spokes are the following edges: 
ei = (vi,«4), e 2 = (v 2l v 6 ), e 3 = (v 3 ,v B ) in (K 5 -l)*, and e; = (^,^+3), « = 1,2,3, 
in K 33 (see Figure [2]). Note that the edges e 2 and e 3 do not intersect in (K 5 — 1)* 
and intersect in K 33 . The described form of K 33 is the graph Q4 of [DG99] . 

Besides the planarity, the graph (K5 — 1)* differs from K 33 by the number of cuts 
of cardinality three. All cuts of cardinality 3 of K 33 are the six one-vertex cuts, 



THE DECOMPOSITION OF THE HYPERMETRIC CONE INTO L-DOMAINS 



15 




Cographic matroid of K 33 Cographic matroid of (K 5 — 1)* 

FIGURE 2. Two cographic matroids 

i.e. cuts containing three edges incident to a vertex. The graph (K5 — 1)*, besides 
the six one-vertex cuts, has a separating cut containing three non-adjacent edges 
one of which is a spoke. (These are the edges (^2,^3), (v5,Vq) and the spoke 
(vi, i> 4 ) of the above description of (K 5 — 1)*.) 

Proposition 4. A maximal by inclusion unimodular subset U of D(4), U 7^ Qi, 
i = 1,2,3, is obtained by a deletion from D{4) of a triple t = (ri,r 2 ,r 3 ) where 
Ti G Qi, z = 1,2,3. Then 

(i) if t has rank 2, then U represents the cographic matroid of the graph K 33 ; 

(ii) ift has rank 3, then U represents the graphic matroid of the graph K 5 — 1, 
which is isomorphic to the cographic matroid of the dual graph (K 5 — 1)*. 

Proof. Since Qi $Z U, for i = 1, 2, 3, U does not contain a triple t. We show that 
a deletion from -D(4) of any triple gives a unimodular set. 

Denote by V(r) the set of all triples of rank 2 containing a root r 6 D(4), If 
r G Qi and r' € Qj, then V(r) fl V(r') is the unique triple of rank 2 containing r 
and r' if i 7^ j, and \V(r) fl V(r')| = if i — j and r 7^ r' . 

Let Ti G Qi, i = 1,2, 3, be the roots of the deleted triple t. Then 

|V(r<)| = 4, |y(ri) n V( rj )\ = 1, and |F(n) n V(r 2 ) f] V(r 3 )\ = or 1. 

Here, whether or 1 stays in the last equality depends on the triple t has rank 3 
or 2, respectively. By the inclusion-exclusion principle, we have 

3 

\ul 1 V(r l )\ = J2\V(r l )\- E l^') n %)l + \r? i=1 V(r i )\. 

i=l l<i<j<3 

Hence, we have the following two cases. 

(i) If I nf =1 V(ri)\ = 1, i.e. if the triple t has rank 2, then | Uf =1 V(n)\ = 
3-4-3-1 + 1 = 10. 

(ii) If I nf =1 V{n)\ = 0, i.e. if the triple t has rank 3, then | Uf =1 V(r t )| = 9. 
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Let V = U re £)(4) V(r) — Uf =1 V(rj). Then V is the set of triples of rank 2 that 
do not contain the deleted roots r 1; r 2 and r 3 . Using the set V we construct a 
cubic directed graph G as follows. We take Ro = D{A) — t to be the set of 9 
arcs of G. As a set of vertices of G we take a subset Vq C V of six triples. An 
arc r G i?o is incident to a vertex v G Vo if the root r belongs to the triple v of 
rank 2. Hence, each vertex u is incident to three edges, and the graph G is cubic. 
Directions of arcs are chosen as follows. Let the roots r,s,p of a triple v satisfy 
the equation r — (r T s)s — (r T p)p = 0. Suppose that the arc r comes in v. Then 
the arc s goes out or comes in the vertex v if r T s = 1 or r T s = — 1, respectively. 
The same assertion is true for the arc p. 

We show that each root r G Ro belongs to exactly two triples of Vq. Let 
r G Q k , then \V(r) PI V(r fc )| = and \V(r) n Vfa) | = |V(r) n VfaOl = 1, where 
{ijk} = {123}. We have two cases: 

(4) V(r) n V(n) = V(r) n Vfa) or V» n V(n) ^ V(r) n Vfa)- 

Suppose that the inequality in (T4J) holds for r. Since | V(r) | = 4, only two 
triples from V(r) belong to V. 

Now, if the equality in (T4J) holds, then it implies that r = (r T rj)rj + (r T rj)rj. 
If t has rank 2, then, up to sign, this gives r = r^, which contradicts to r G Ro- 
Hence, if t has rank 2, the inequality of (jl]) holds for all r E R$. We can set 
Vo = V, since V contains 16 — 10 = 6 triples. It is easy to verify that G is 
isomorphic to K 33 . 

If t has rank 3, then r = r^, where is defined in equation ([3]). The roots r^, 
ij = 12, 23, 31 when the equality holds in|4j form a triple vq of rank 2. Obviously, 
■uo G V. We set Vo — V — {vo}. Since the root belongs to two triples vq and 
(rij,ri,rj), the remaining triples of V(rij) belong to Vo. We saw that each root 
r G Ro, r 7^ for ij = 12,23,31, belongs to two triples of i?o- Hence, in the 
case when t has rank 3, the graph G is well defined. It is easy to verify that G is 
isomorphic to (K 5 — 1)* and v corresponds to a separating cut of cardinality 3. 

Arcs of G are labeled naturally by roots from the set R^. This labeling gives a 
representation of the cographic matroid of G by vectors of i? - 

Note that triples of rank 2 are equivalent under action of the automorphism 
group of Py(D 4 ). Similarly, all triples of rank 3 are equivalent under the au- 
tomorphism group of Py(D4) extended by changing signs of roots. Hence, any 
explicit representations of the cographic matroids of the graphs K33 and (K5 — 1)* 
for fixed triples of rank 2 and 3 proves that above labeling give representations 
for all pairs of triples of rank 2 and 3. □ 

6.4. L-types in £(HYP 5 ). From Proposition H] we deduce that each L-domain 
corresponding to D4 is contained in 64 = 4 3 different L-types. 48 of them are of 
type II and 16 are of type III. 

Recall that PviD^) is free along lines spanned by roots of the root system 04^, 
vectors of which are parallel to diagonals of the cross-polytopes f3& of the Delaunay 
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partition of the lattice D 4 . Note that bases related to the forms a(D 4 ) £ £(HYP 5 ) 
contain a diagonal of a /3 4 . 

Lemma 5. Let the basis related to a(D 4 ) contains a diagonal q £ D 42 of a 
cross-polytope /3 4 . Then the L-domain of the parallelohedron P q = Py(D 4 ) -f z{q) 
related to this basis, i.e. the L-domain of the form a(D 4 ) + \ f q , does not belong 
toC(HYP 5 ). 

Proof. The parallelohedron P q has a non-zero width along the line l(q) parallel 
to q. The lattice L q of P q has a lamina H which is orthogonal to q. The lamina 
H separate the cross-polytope /3 4 with a diagonal q into two Delaunay polytopes, 
each being a pyramid with a base (3% orthogonal to q and lying in H. These two 
pyramids have the end-points of the diagonal q as apexes. Hence vertices of the 
basic simplex of a(D 4 ) + \ f q belong to two distinct Delaunay polytopes. This 
implies that the L-domain of L q does not belong to ^(HYP 5 ). □ 

Proposition 5. The cone ^(HYPs) contains 12 principal L-domains, 120 L- 
domains of type II and 40 L-domains of type III, total 172 L-domains. 

Proof. By Proposition [2J £(HYP 4+1 ) contains |4! = 12 principal L-domains. 

The closure of an L-domain of types II and III is the convex hull of a(D 4 ) and 
a dicing facet F{U) related to a unimodular subset U C D 42 of 9 vectors that 
are free for Py(D 4 ). 

Each tile of the Delaunay tiling of the root lattice D 4 is the regular four- 
dimensional cross polytope /3 4 . Any affine base of /9 4 contains exactly the two 
vertices w, w' of a diagonal and 3 vertices t>i, i> 2 , ^3 chosen from the other diago- 
nals. We have dp 4 (w,w') = 4 and 2 for all other pairs. There are = 10 ways 
to choose a pair {w,w'} in a 5 elements set so there are exactly 10 rays «j(D 4 ) 
representing D 4 in HYP 5 . 

Every L-domain aj(D 4 ) is contained in the closure of 64 primitive L-domains 
T>(U), but not all of them are included in ^(HYPs). Each L-domain T>(U) has 
the form T>(U) = com;(aj(D 4 ) + F(U)), where the subset U C D 42 is obtained 
by a deletion of a triple from D 42 . By Lemma the inclusion T>(U) C ^(HYP 5 ) 
implies U does not contain the diagonal q which is a basic vector of aj(D 4 ). Let 
q — 7*1 £ Qi, then we have 4 choices for the root r 2 in Q 2 . For the third root either 
we have a rank 2 triple and r% is completely determined or we have a rank 3 triple 
in which case there are 3 choices. This means that a ray a.;(D 4 ) is contained in 4 
L-domains of type II included in ^(HYPs) and 12 L-domains of type III included 
in £(HYP5). This means that there are 120, respectively 40 L-domains of type 
II, respectively III in £(HYP 5 ). □ 
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